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ABSTRACT 

We study the evolution of galactic bars and the link with disk and spheroid formation in a sample of 
zoom-in cosmological simulations. Our simulation sample focuses on galaxies with present-day stellar 
masses in the 10 10-11 M Q range, in field and loose group environments, with a broad variety of mass 
growth histories. In our models, bars are almost absent from the progenitors of present-day spirals 
at z > 1.5, and they remain rare and generally too weak to be observable down to z « 1. After this 
characteristic epoch, the fractions of observable and strong bars raise rapidly, bars being present in 
80% of spiral galaxies and easily observable in two thirds of these at z < 0.5. This is quantitatively 
consistent with the redshift evolution of the observed bar fraction, although the latter is presently 
known up to z w 0.8 because of band-shifting and resolution effects. Our models hence predict that 
the decrease in the bar fraction with increasing redshift should continue with a fraction of observable 
bars not larger than 10-15% in disk galaxies at z > 1. Our models also predict later bar formation 
in lower-mass galaxies, in agreement with existing data. We find that the characteristic epoch of 
bar formation, namely redshift z w 0.8 — 1 in the studied mass range, corresponds to the epoch at 
which today's spirals acquire their disk-dominated morphology. At higher redshift, disks tend to be 
rapidly destroyed by mergers and gravitational instabilities and rarely develop significant bars. We 
hence suggest that the bar formation epoch corresponds to the transition between an early "violent" 
phase of spiral galaxy formation at z > 1 and a late "secular" phase at z < 0.8. In the secular 
phase, the presence of bars substantially contributes to the growth of the (pseudo-)bulge, but the 
bulge mass budget remains statistically dominated by the contribution of mergers, interactions and 
disk instabilities at high redshift. Early bars at z > 1 are often short-lived, while most of the bars 
formed at z < 1 persist down to z = 0, late cosmological gas infall being necessary to maintain some 
of them. 

Subject headings: Galaxies: formation — Galaxies: evolution — Galaxies: spiral — Galaxies: structure 
- Galaxies: bulges 
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1. INTRODUCTION 

Bars are one of the most frequently and easily quanti- 
fied substructures in spiral galaxies, and are hence often 
used as a tracer of galaxy evolution. Most spiral galaxies 
today contain a central bar, although with largely var i- 
able amplitudes (jBlock et all I2002t IWhvte et al.ll200l . 
Spiral arms are equally common in optical light, but are 
much weaker in the near-infrared light that more closely 
traces the stellar mass distribution, and the strength of 
bars is generally easier to quantify independently of the 
imaging sensitivity. 

Both the formation of bars and their time evolution 
are connected to the baryonic and dark matter properties 
of their host galaxies, and their mass assembly history. 
Bars form spontaneously in stellar disks that are suffi- 
ciently massive and dynamically cold to be gravitation- 
ally unstable, with typical Toomre stability parameters 
Q k 1.5 - 2.0 (lToomrd[T96l ICombes fc Sandersl HMD: 
iGerin et all 119901 : iCombes fc Elmegreenlll993f r and can 
also be amplified by dynamical friction on the dark mat- 



ter halo (| Athanassoulal 120021 ) . The gaseous content can 
also trigger bars: gas helps to form outer spiral arms, 
which can remove angular mome ntum from the inner re - 
gions and strengthen a bar seed (|Bournaud et al.|[2005n . 

Once formed, bars evolve through the exchange 
of angular momentum with the dark mat ter halos 
(jWeinbergl [19851 IDebattista fc Sellwoodl 1200(1 . as well 
as with stellar and gaseou s disks dFriedli fc Benzj [l99l 
IBournaud fe Combed 12001 IBournaud et al.ll2005f) . This 
can reinforce bars when they lose angular momentum, 
but c an also weaken and d estroy them in the opposite 
case (IBournaud et al.ll2005l ). B ars can also be weakened 
or des troyed by the growth of central mass concentra- 
tions (Hasan fe Normal 1 19901: IPfenniger fe Friedhlll991t 
IHasan et al.lll993D . however central concentrations with 
low-enough mass and/or low-enoug h mass growth rates 
could have little effect on real bars (|Athanassoula et al.l 
2005). If bars are formed in conditions where they 
are intrinsically short-lived, sufficient accretion of ex- 
ternal gas onto the disk could enable their survival 
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or re- formation (iBournaud fc Combes! [2002). In ad- 
dition, galaxy interactions could in theory trigger bar 
(re-)formation (iGerin et al.lll990HMiwa fc Nogu"chlll998l: 
iBerentzen et alJ 120041 ) . although observations do not 
show a clear environmental dependence of the bar in 
disk galaxies (Ivan den Berghl|2002t lAguerri et al.l [20091 : 
iBarwav et al1l201lHLee'eT^[20T2[ ). 

Hence, the fraction of barred galaxies, and the red- 
shift evolution of this fraction, are fundamental tracers 
of the evolution history of galaxies: this indicates when 
disks became sufficiently massive and self-gravitating to 
be bar-unstable, and whether the conditions for bars be- 
ing long-lived or re-formed were met. Furthermore, bars 
can directly drive structural evolution of their host galax- 
ies. They trigger radial gas flows and m ay provide gas 



to nuclear disks and ce ntral black holes ( Knayen et all 
120001 lLaine et al.ll2002t lLaurikainen et al.l 120041 ) . Thev 



can also thicken through vertical resonances, leading to 
the formation of pseudo-bulges - i.e. , bulges with rela- 
tively low concentration and substantial residual rotation 
(e.g.. Bureau fe Freemanl 119991: iKormendv fe Kennicuttl 



12004 [Martinez- Valpucsta!eLaij[2C^|) . 



In the nearby Universe, the bar fraction in disk galax- 
ies is very high. Depending on classification tech- 
niques, the fracti on of strong bars in the optical light 
is at least 50% (jBarazza et alJ |2008[) . Optical clas- 
sifications reveal roughly one third of strongly barred 
galaxies, one third of weakly or moderately barred 
gala xies, and one third of op tically unbarred galax- 
ies (jde Vaucouleurs et alJ 119911 ) . In the near-infrared, 
where weak bars are not obscured by dust and more 
easily distinguis hed from spiral arm s , the bar fraction 
is at least 80% (lEskridge et al.ll200dt Block et alJlgOOl 
iWhvte et aLll2002t IMenendez-Delmestre et al.H2007t) . 

The first searches for bars at redshift z > 0-5 fo und a 
very low bar fraction (e.g., lAbraham et alJ fl999l ) . pos- 
sibly because of small number statistics. Their work 
also illustrated the difficulties to identify bars at high 
redshift: the observed optical light traces the ultraviolet 
emission, in which bars are harder to detect, even locally. 
Near-infrared data revealed a number of barred galax- 
ies at z > 0.7 (ISheth et~al~ll2003l lElmegreen et"aLll2004l : 
Uogee et al.ll2004|) . The first sample large enough to ro- 
bustly quantify the redshift evolution of the bar fraction 
without being affected by resolutio n and band-shifting 
bias up to z w 0.8, was studied by ISheth et alJ (j2008f ) 
in the COSMOS field. These observations indicate that 
the bar fraction drops by a factor of about three from 
z = to z — 0.8. This result holds both for all observ- 
able bars and for strong bars separately, and using cither 
visual cla ssifications or quan titative estimates of the bar 
strength. ISheth et alJ (|2008f ) also found a downsizing-like 
behaviour for bar formation, i.e. more massive galaxies 
tend to get barred at higher redshifts . This trend can ex - 
plain why previous studies, such as Uogee et alJ (|2004D . 
using shallower data targeted to more massive systems, 
observed higher bar fractions - but still consistent wit h 
a declining bar fraction (see also lElmegreen et al.ll2004l ). 

In this paper, we study the evolution of bars in a sam- 
ple of cosmological zoom-in simulations of 33 galaxies 
with present-day stellar masses ranging from 1 x 10 10 to 
2 x 10 11 M.q, in field and loose group environments. The 
simulation technique and structural evolution of these 



galaxies (bulge and disk fra ctions, angula r mom entum 
evolution) were presented in IMartig et all (|2012l ). The 
paper is organized as follows. In Section [2] we present 
our simulations and methods for the identification of bars 
and morphological analysis. In Section [3] we analyze the 
redshift evolution of the bar fraction in the whole sam- 
ple and in disk-dominated galaxief], using quantitative 
measurements of the bar strength. Our main result, the 
emergence of bars along the cosmic time that traces the 
epoch of thin disk formation subsequent to the growth of 
spheroids and thick disks, is presented in Section [H Sec- 
tion [5] studies the lifetime of bars and its dependence on 
external gas accretion. In Section [51 we quantify the con- 
tribution of bars in the late growth of bulges, comparing 
to unbarred galaxies. Finally, we discuss and summarize 
our results in Sections [7] and [8j respectively. 

2. SIMULATIONS AND ANALYSIS 

2.1. Simulation Sample 

The simula t ion sa mple studied here was presented in 
IMartig et all (pOll) . It comprises 33 field and loose 
group galaxies modeled at 150 pc resolution from redshift 
5 down to redshift zero with present-day stellar masses 
ranging from 1 x 10 10 to 2 x 10 11 M Q . 

The zoom-in s imula t ion t echnique is fully described 
in IMartig et all (|2009l 120121 ). Dark matter haloes are 
selected in a large volume, dark matter-only simula- 
tion, performed using A-cold dark matter cosmol ogy with 
the A daptive Mesh Refinement code RAMSES (|Tevssierl 
[2001 . The zoom-in simulation for each selected halo is 
performed by recording and replicating the mass inflow 
through a spherical boundary at the final virial radius of 
the selected halo. The boundary conditions used for the 
zoom-in simulation hence replicate all minor and major 
mergers as well as diffuse infall, as imposed by the initial 
cosmological simulation. 

Our zoom-in simulations use a spatial resolution of 
150 pc, and a mass resolution of 10 4 ~ 5 M Q (1.5 x 10 4 M Q 
for gas, 7.5 x 10 4 M Q for stars and 3 x 10 5 M for dark 
matter particles). The code used to model gravity for 
gas, stars a nd dark matter is th e part i cle-me sh code 
described in IBournaud fe Combesl ([20021 120031 ). Inter- 
stellar gas dynamics is modeled with a sticky-particle 
scheme (studied in the context of bar evolution in 
IBournaud et al.l I2005D , which has the drawback of ne- 
glecting thermal pressure (especially in hot halos, not 
crucial in the studied mass range), but the advantage of 
modeling the properties of supersonic turbulent pressure 
in cold gas phases - which is the physically dominant 
pres sure term in t he star-forming interstellar medium 
(e.g. iBurkertl 120061) and is hardly modeled by other hy- 
drodynamic techniques unless resolutions of 1 — 10 pc and 
low thermal cooling floors are reached dTevssier et all 
[2011 IBournaud et^[201lT: iHopkins et al.l 120 111 ). Star 
formation is com puted with a Schmidt-Kcnnicutt law 
(|Kennicuttl I1998D with an exponent of 1.5 and an ef- 
ficiency of 2%. The star formation threshold is set 
at 0.03 M Q pc~ 3 . Energy feedback from supernovae ex- 
plosions using a kinetic scheme as well as the contin- 
uous gas mass-loss from stars (jJungwiert et alJ 120011 : 
IMartig fc Bournaudl l2"Ol"of ) are included. The z = 5 

1 most galaxies in this sample are disk-dominated spirals at z = 
0, but a larger fraction is spheroid-dominated at z > 1 
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seed for the central galaxy and the incoming compan- 
ions are implemented with arbitrary disk+bulgc models; 
tests have shown that the seed properties at z = 5 have 
no substantial impact on the structural evolution from 
z = 2 out to z = 0, because rapid evolution and mass 
growth at z > 3 washes out the initial assumptions. 

In our study, we focus on the redshift range from 
z = 2 to z = 0. Detailed descr i ption of analyzed sam- 
ple is presented in IMartig et al.1 ([2012D . Here, we recall 
typical masses and star formation rates (SFR) at dif- 
ferent redshifts. At z = 0, typical mass and SFR of 
simulated galaxies are 4 — 5 x 1O 1O M0 and 3M yr _1 , 
at z = 1, 2 x 10 10 M Q and 5M Q yr~ Y and at z = 2, 
2 — 3 x 1O 9 M and 10M Q yr _1 , respectively. Hence, a 
broad rang e of mass growth h istories is covered by the 
sample (see IMartig et al.ll2012l for details). 

2.2. Bar analysis 

The method used to determine the presence, length 
and strength of a bar is based on the azimuthal spectral 
analysis of surface density profiles. This is obtained by 
considering the stellar surface density of each galaxy in 
polar coordinates, decomposed into Fourier components 
in the form: 

S(r, 9) = E (r) + A ™M cos ( md ~ $ ™( r ))> (1) 

m 

where E(r, 6) is the stellar surface density, 8 the az- 
imuthal angle given in rotating frame in which the bar is 
fixed and r the radial distance. A m and $ m are the asso- 
ciated Fourier amplitude and phase, respectively. Eo(r) 
is the azimuthally-averaged profile of the stellar surface 
density. The analysis is performed on a face-on projec- 
tion: the spin axis of the whole stellar content of the 
galaxy is used to define the corresponding line of sight. 
The center of mass of the stars within the central 10 kpc 
is taken to be the center of the galaxy for the Fourier 
analysis. 

A typical signature of the presence of bar is the promi- 
nence of even components, especially m = 2, within the 
bar region. The identification of bar is possible by study- 
ing the phase <&2(r) which is constant with radius in the 
bar region, as opposed to a two-armed spiral mode ($ m 
varies linearly with radius for an m- armed spiral mode). 

After careful examination of the whole sample at dif- 
ferent redshifts, we decide that a bar is present if §2(1") is 
constant to within ±5 around the median value over a 
large-enough region, hereafter the "bar region" . This bar 
region should start at radii between 900 pc and 1500 pc, 
and cover a radial range of at least 1500 pc. We search 
the starting point of the bar region at r > 900 pc be- 
cause central asymmetries or off-centering in the central 
resolution elements can cause variations of $2 at smaller 
radii even for visually barred systems, while starting the 
search at 900 pc is found to never exclude systems that 
are visually identified as barred. We stop the search at 
1500 pc as no bar identified visually starts its $2 — con- 
stant region at larger radii. We impose the bar region 
($2 constant to within 10°) to cover a radial range of 
at least 1500 pc because this selects all bars identified 
visually and this excludes spiral arms, as the latter typ- 
ically have a variation of $2 of a few tens of degrees 
over a few kpc. The choice of cutoff in the bar length 
is given mainly by the resolution. We require the bar to 



lie within at least three resolution bins used in Fourier 
analysis. Moreover, typical lengths of sign ificant bars are 
in general > 2 kpc ((Barazza et all l200cf) . Smaller bars 
are usually either nuclear bars (which are not the focus 
of this study) or weak ones. We could possibly miss some 
of these short and weak bars, but this does not have any 
impact on our conclusions on bars that would be strong 
enough to be observed at z > 0. 

Once a bar region meeting to the previous criteria is 
found, the galaxy is classified as barred and the bar 
strength S is m easured following the definition proposed, 
for instance, bv lAguerri et al.1 (|1998f ): 



S 



= rr 1 

bar 



-A 2 
-j-dr, 
A o 



(2) 



where ri> ar is the outer radius of the bar region. 

Figure [1] shows two examples of barred and unbarred 
galaxies and the corresponding radial profiles of $2- 
With this technique, the automated identification of 
barred galaxies is in agreement with all the cases that 
have been visually examined. We classify bars into strong 
bars and weak bars by considering the strength S (Eq. [5]) . 
In the following, we will distinguish: 

• all detected bars, including very weak and short 
bars, 

• observable bars, with a strength S > 0.2, which is 
the typical detection l imit used in obser vations up 
to high redshift (e.g., iSheth et all 120081 ) - even at 
z = 0, weaker bars may be confused with spiral 
arms unless the $2 phase can be probed accurat ely 
with very deep imaging (e.g., IBlock et all 12002^ . 

• and strong bars with S > 0.3 unless specified oth- 
erwise. 

At this stage, the identification of bars can still be con- 
fused with flattened early type galaxies, especially for 
weak bars. A spheroid-dominated galaxy, in the face- 
on projection defined by the stellar spin axis, can have 
flattened isophotes in its central regions, which could be 
identified as a "bar" with our definition, while this corre- 
sponds to a triaxial (part of) stellar spheroid rather than 
a real bar: such an example is shown on Figure [2] This 
occurs only with weak bars, and there are few spheroid- 
dominated galaxies in our sample at z = 0, so the prob- 
lem is unimportant there, but such cases are somewhat 
more frequent at increasing redshift. We thus need to dis- 
tinguish such "fake" bars from real (if weak) bars. We 
find that examining two edge-on projections of the stel- 
lar mass density unambiguously solves the problem for 
our whole sample. Real weak bars are found in galaxies 
with a massive disk component, hence edge-on projec- 
tions of the stellar density are substantially flat. "Fake" 
bars are found in galaxies that are spheroid-dominated 
(at least in the central few kpc) and the edge-on projec- 
tions are quite round (see Figure [3] for an illustration). 
We thus decide that galaxies for which the strength S of 
the m = 2 mode in two orthogonal edge-on projection^ 
is greater or equal to 0.3 correspond to real weak bars in 

2 In these projections the m = 2 mode strength is not a bar 
strength, but traces the presence of an edge-on disky component 
in the stellar mass distribution. 
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Fig. 1. — Examples of the fitting method results (left panels) and the corresponding surface density maps (right panels) viewed face-on 
(50 X 50 kpc 2 ) for two galaxies in our sample. The top panel shows a (strongly) barred galaxy with the region where the phase $2 is 
almost constant identified as the bar region, while the bottom panel corresponds to an unbarred galaxy for which no region of constant <3?2 



is detected. The color coding scale of surface density map is logarithmic going from 
dark red. 
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Fig. 2. — Example of a possible "fake" bar. A region of constant $2 is identified at r < 2.7 kpc (left panel), corresponding to the detection 
of a short bar. However, the face-on density map (right panel, 50 X 50 kpc 2 ) suggests that the identified region could instead be the core 
of an elliptical, i.e., a "fake" bar, which is confirmed by the examination of edge-on projections performed in Figurc[3] The color coding of 
the projected density map is the same as in Figure [T] 
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Fig. 3. — Examples of surface density maps (50x50 kpc 2 ) in three different projections for two simulated galaxies. The corresponding 
rcdshift is z = and 2 = 1 for galaxies in top and bottom panels, respectively. If a potential bar is identified in the face-on projection, the 
two orthogonal edge-on projections are used to discriminate real bars from triaxial cores in spheroids or "fake" bars. The color coding of 
the projected density maps is the same as in Figure [T] 



disky stellar systems. Other cases are almost round in all 
projections (S < 0.3 in all projections) and are consid- 
ered to be spheroid-dominated galaxies with "fake" bars, 
i.e. moderately flattened central regions. 

These "fake" bars are rare at low redshift (about 6% of 
all galaxies) and represent less than 15% of all galaxies 
at z > 1.5. Visual inspection of the results showed that 
at most two thirds of the "fake" bars are successfully re- 
moved by the technique described above at any redshift, 
so any remaining contamination of the bar fraction would 
be quite minor. Moreover, most of the results studied 
in this paper relate to the fraction of bars among disk- 
dominated galaxies at various red shifts (based on a Sersic 
index measurement (|Sersidll963() . see Section [273]) . which 
is not contaminated by "fake" bars. Indeed, these "fake" 
bars are flattened spheroid-dominated galaxies, which 
are naturally removed from the Sersic index-selected 
sample of disk-dominat ed galaxies. 

It has been shown by iBournaud et all (|2005l ) that the 
m = 2 modes are not strongly affected by the dissipa- 
tion parameters used in sticky particle codes. We have 
performed additional checks to estimate the impact of 
the centering of the galaxy on the result of the Fourier 
analysis. Changing the in-plane coordinates of the cen- 
ter used for the Fourier decomposition by as much as 
300 pc, wc do not find any significant change in the de- 
tection and strength of bars classified as observable (as 
defined above). Some galaxies identified as unbarred can 
appear weakly barred, and weak bars may not be de- 
tected anymore, when the center used for the Fourier 
decomposition is moved by as much as 300 pc. Such am- 
biguous cases are quite rare for centering offsets smaller 
than 300 pc. The overall rate of misclassification due to 



the galaxy off-centering is not larger than 15% and af- 
fects mostly bars that would not be observable at high 
redshift. 

2.3. Morphology analysis: Identifying disk-dominated 
and bulge-dominated galaxies 

In order to distinguish disk-dominated galaxies from 
earlier- type ones, we fit the radial profile of the stellar 
mass surface density with a Sersic profile of form: 



£(r) = ao exp 



r 

-ai \ — 

Jo 



(3) 



where ao, a\ are normalization constants, n is the Sersic 
index, r is the radius and t$ is the scale-length. Galax- 
ies with a Sersic index n < 2 are classified as disk- 
dominated. The fitted range is r 50 < r < 2 x r 90 for 
unbarred galaxies and bar < r < 2 x rgo (or bar < r < 
2.5 x i 2 5 if the length of the bar is shorter than rgo) 
for barred galaxies, where r$o and rgo are radii contain- 
ing 50% and 90% of mass, respectively, bar is the cor- 
rected length of the bar (length of the bar - 900 pc) and 
i25 the 25 th isophote. The scale-length ro is set to the 
value of rso for each galaxy. The ranges for our set of 
simulated galaxies are chosen to produce a satisfactory 
match to the observed morphology of these galaxies on 
one side and staying relatively simple with comprehensi- 
ble physical interpretation on the other side. Indeed, the 
present classification is reproducible in observations and 
it agrees with the classification based on a decomposition 
of disk and bulge co mponents in the 6-D phase space by 
iMartig et all l)2012l ) for the vast majority of the sample. 
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The galaxies presently classified as disk-dominated have 
bulge-to-total mass fraction below 35%. 

Note that the identification of a disk component to re- 
ject "fake bars" (Section I2.2[) was based on a different 
criterion. There we aimed at probing the presence of a 
substantial disky component hosting the bar, without re- 
quiring it to dominate the stellar mass distribution, and 
we used edge-on projections for this. Here we use the 
Sersic index of the face-on projection to probe whether 
the whole galaxy is dominated by a disk component typ- 
ically twice more massive (at least) than the bulge. As 
the Sersic profile is fitted to the stellar mass density, it 
is not sensitive to the nuclear concentrations of star for- 
mation. Nevertheless, we checked the consistency of the 
identification with the bulge-to-total ratios obtained by 
iMartig et all (poll ). 

Figure |4] shows the evolution with redshift of the frac- 
tion of disk-dominated galaxies within all 33 simulated 
galaxies. 70% of galaxies are found to be disk-dominated 
at z = 0, and th eir bulge fraction is 23% (from the 
IMartig et al.l 120121 decomposition), consistent with clas- 
sifying them as "spiral galaxies" . The remaining 30% are 
galaxies with an average bulge-to-total ratio of 0.43 and 
could be classified as early-type SO galaxies. We find a 
decreasing fraction of disk-dominated galaxies as a func- 
tion of redshift, down to a fraction of 20% at z = 1. 
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Fig. 4. — Evolution with redshift of the fraction of galaxies clas- 
sified as spirals according to their Sersic index as defined in the 
text. At z = about two thirds of galaxies are found to be disk- 
dominated while for z > 1 this fraction is ~ 0.2 with ±1<t fluctua- 
tions. Here and in all following figures, the error bars on measured 
fractions are Poissonian, the average error bar being attributed to 
zero values, unless stated differently. 

3. THE REDSHIFT EVOLUTION OF THE BAR FRACTION 

3.1. Bar fraction history 

Figure[5]shows the redshift evolution of the bar fraction 
among galaxies classified as spirals (i.e., disk-dominated). 
The top panel shows this fraction including all values of 
the bar strength, while the middle and bottom panels 
show the redshift evolution of the bar fraction for observ- 
able bars (strength S > 0.2) and strong bars (S > 0.3), 
respectively. 

The total bar fraction among spiral galaxies declines 
with increasing redshift. At redshift zero, 80-90% of spi- 



ral galaxies contain a bar, while at z ~ 1 this fraction 
drops to about 50%, and to almost zero at z ~ 2. Sim- 
ilarly, the observable bar fraction and strong bar frac- 
tion decline from about 70% at z — to 10 — 20% at 
z ~ 1. Observable and strong bars are virtually absent 
at z > 1.5. 

The steady decline of the total, observable and strong 
bar fractions from z = 0toz^l — 2 cannot be attributed 
to a decline of the fraction of disk-dominated galaxies, 
since the bar fractions mentioned above are measured 
among disk-dominated (spiral) galaxies. There is a de- 
cline in the fraction of disk galaxies with increasing red- 
shift (FigurcHJ), and as a result the bar fraction among all 
galaxies declines even more rapidly with redshift (from 
«60% at redshift z — to less than 5% at z = 1, Fig- 
ure [6]) but the decline of the bar fraction occurs even 
considering only disk-dominated galaxies. 

We compare in Figures [7] and |8] our data to the ob- 
servations of the z = — 0.8 sa mple in COSMOS and 
the SDSS by iSheth et all (|2008f ). which is the largest 
and highest-redshift dataset measuring the observed bar 
fraction to date. The redshift evolution of the bar frac- 
tion is in quantitative agreement between models and 
data, considering both the fra ction of obs e rvable bars 
and of strong barfl Note that ISheth et all (|2008f ) used 
an ellipse-fitting technique to measure bar strengths, 
but such techniques give quantitative results quite con- 
sistent with Fourier de compositions (jBlock et al.l 120021 : 
lLaurikainen et aLll2002f) . as well as with visual estimates 
bv ISheth et al.1 (|2008| ) themselves. We also note that the 
fraction of bars that could be robustly observed at z > 0, 
i.e. with a strength larger than 0.2, is lower than the total 
fraction of bars identified in the simulations including the 
weakest bars, especially at redshift z ~ 1 (see top panel 
of Figure [5] versus Figure [7]) , but both fractions follow a 
similar evolution with redshift. 

Observations of the bar fraction are typically limited 
to re dshift z < 0-8 to date (with very rare z ~ 1 exam- 
ples, |Ehnegreen~£tjd- 2004). Our models agree with the 
observed trends up to those redshifts, and suggest that 
the observed decline should continue at higher redshifts, 
observable bars being almost absent at z > 1. 

3.2. Dependence on stellar mass 

Figure shows the redshift evolution of the total bar 
fraction among spiral galaxies for low-mass and high- 
mass systems separately. We consider galaxies above and 
below the median mass of the sample in each redshift bin. 
The bin size was adjusted so that the Poissonian error in 
each bin does not exceed 20%. 

We find that bars form later in lower-mass galaxies, 
in a "downsizing" -l ike mode, which i s consistent with 
the observations of ISheth et al.l (|2008f ). These authors 
proposed that merging activity, which is more common 
at high redshift, could affect low- mass systems more 
severely by heating them and thus delaying or prevent- 
ing bar formation. However, mergers alone cannot ac- 
count for the observed trend in our models. First, ma- 
jor mergers of spiral galaxies often result in the forma- 
tion of spheroid-dominated galaxies, even at high red- 
shift (jBournaud et al.ll2011h . We find that the majority 

3 defined in Figure \E\ as bars st ronger than 0.4 (S > 0.4) for 
consistency with ISheth et al.l l|2008fl 
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Fig. 5. — Evolution with redshift of the total bar fraction in- 
cluding even weak bars (top panel), observable bar fraction (mid- 
dle panel), and strong bar fraction (bottom panel) among spiral 
galaxies (i.e., galaxies with a Sersic index n < 2 in each redshift 
bin). 

of disk galaxies are barred today, and bar formation is 
found in our sample even when considering only disk- 
dominated galaxies with low Sersic indices, so the pro- 
cess should be independent on the occurrence of major 
mergers. Second, minor mergers and interactions can 
cause the destruction or formation of a bar, depending in 
particular on the orb i tal pa rameters (|Gerin et al.l 119901 
iBerentzen et~aT1l20"03l [200l . Since bars at high redshift 
(z > 1) are rare, any net effect of mergers/interactions 
should be bar formation in previously un-barred systems, 
rather than bar destruction in previously barred systems. 

Over the entire sample, our results suggest that the 
epoch of bar formation is the typical epoch at which 
galaxies start to be dominated by a kincmatically cold, 
thin stellar disk (see Section [4}. A possible explanation 
of the "downsizing" of bar formation could be that these 
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Fig. 6. — Evolution with redshift of the fraction of barred galaxies 
among all galaxies, where only strong bars (S > 0.3) are consid- 
ered. Note that the bottom panel of Figure |51 showed strong bars 
among galaxies classified as disk dominated (Sersic index n < 2). 
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Fig. 7. — Comparison of the redshift evolution of t he bar frac- 
tion i n spiral galaxies with COSMOS and SDSS data llShcth et al.l 
[2001) . The lower limit on the bar strength of simulated galaxies 
is set to 0.2, so that onl y observable bars a re shown. The error 
bars are calculated as in [Shcth et al. (2008), assuming binomial 
statistics. 



modern spiral disks themselves form later in lower-mass 
galaxies. The continuation of rapid mass accretion onto 
lower-mass systems down to lower redshift could keep 
their disk violently unstable, with giant clumps and ir- 
regular structures rather than bars, as further discussed 
in Section [7] More massive spiral galaxies seem to be 
large ly in place and alre ady disk-dominated at z ~ 1 
(e.g.. (Sargent et al.ll2007l) . 

We also find that at intermediate redshift (z < 0.6) the 
fraction of bars in low- mass galaxies continues to increase 
with decreasing z, while in high-mass systems, it stays 
roughly constant from z ~ 0.6 to z = at a relatively 
high value of ~ 70%. At high redshift (z > 1), the most 
massive disk galaxies do contain some bars (although rel- 
atively rare, and in general too weak to be easily observ- 
able), while lower mass disk galaxies are systematically 
unbarred. 
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Fig. 8. — Comparison of the evolution with redshift of the frac- 
tion of strong bars (here considering strengt hs greater than 0.4 ) 
in spirals of COSMOS and SDSS data fro m Shcth cFaTl l2008Tl . 
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Fig. 9. — Evolution with the redshift of the fraction of barred 
galaxies among spiral galaxies according to their mass. In each 
redshift bin, galaxies are classified as having a stellar mass above or 
below the median value, which is 4 X 10 9 M Q at z=2, 2.5 X 10 10 M Q 
at z ~ 1 and 6 x 10 10 M at z=0. 

4. WHEN DO BARS EMERGE? THE EMERGENCE OF 
MASSIVE THIN DISKS 

While the total bar fraction and strong bar fraction 
are both high at low redshift, they become very small or 
close to zero at z > 1. Most of the bars found at z > 1 
in the simulations are quite weak and would not even 
be considered with a bar strength threshold at S = 0.2, 
hence they would be hardly observable at high redshift 
(see Figure [7]) . 

The redshift evolution of the fraction of disk- 
dominatcd3 galaxies in our simulation sample was shown 
on Figure 01 Surprisingly, it is quite similar to the red- 
shift evolution of the fraction of strong bars or moderate 
(observable) bars (see Figure 0] compared to the bottom 
panel of Figure [5] and Figure [7]). This suggests that the 
epoch of bar formation is also the epoch at which galaxies 

4 according to criterion defined in Section 12.31 
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Fig. 10. — Evolution with the redshift of the fraction of spiral 
galaxies among all galaxies according to their mass. The binning 
is same as in Figure [9] 

that are presently spirals start to be dominated by stel- 
lar disks. Indeed, at higher redshift, the progenitors of 
these present-day spirals are sometimes disk-dominated 
but are often spheroid-dominated, or have irregular mor- 
phologies such as clumpy disks (i.e., short lived unsta- 
ble disks), or are interacting and merging systems (see 
iMartig et al.ll2012l for a thorough study of the structural 
evolution of the simulated g alaxies in the s a me sa mple) . 

Our study and that in IMartig et all (|2012f ) indi- 
cate that these present-day spiral galaxies grow mostly 
through two phases. During an early "violent" phase at 
z > 1, the system is quite often disturbed by impor- 
tant mergers (major mergers or multiple minor merg- 
ers) as well as violent disk instabilities (giant clumps). 
The morphology can evolve during this phase from disk- 
dominated to a spheroid without being stabilized toward 
the hnal disk-dominated structure. At z < 1, present-day 
spirals evolve mostly through a "secular" phase when the 
morpholo gy is generally stab ilized to a disk-dominated 
structure ( Martig et al.[|2012l ). The bulge growth is then 
only slow and limited as important mergers and violent 
disk instabilities become rare, and almost absent after 
z ~ 0.7. 

ioii), 



As measured in IMartig et al.l (|2012l ) , there is no cor- 
relation between the morphologies and the disk/bulge 
fractions between the early violent phase (z > 1) and 
the present day. While the systems enter the "secular" 
phase after z w 1 the bulge/disk fractions become more 
and more tightly correlated with the final z = val- 
ues. Stars formed or accreted in the early violent phase 
end-up mostly in the thick disk, stellar halo and bulge 
at redshift zero. Stars formed in the secular phase at 
z < 1 mostly grow the modern thin spiral disk, with a 
substantial but non dominant contribution to late bulge 
growth, especially in barred systems (see Section |6] for 
late secular bulge growth). 

The epoch of bar formation covers a relatively nar- 
row range z ~ 0.7 — 1 for the vast majority of galaxies 
in our sample, although the sample spans stellar masses 
of 1 x 10 10 to 2 x 10 11 Mq. The comparison with the 
structural evolution shows that the (observable) bar frac- 
tion is constantly very low in the early violent phase 
of spiral galaxy formation, a phase during which disk- 
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dominated systems can be present but are generally de- 
stroyed/reformed over short timcscales, and uncorrelated 
with the final disk fraction. Once spiral galaxies en- 
ter their "secular phase" at z < 1, bars rapidly form 
and the bar fraction rapidly raises. This is the epoch at 
which the final thin spiral disk starts to form and domi- 
nate the stellar structure of these galaxies: the disk and 
bulge fractions can continue to evolve down to z = 
but in this phase the disk is not destroyed/reformed 
anymore, although the bar itself may sometimes be de- 
stroyed/reformed. At this point, the formation of a long- 
lived massive thin disk allows the bar to form in most of 
the progenitors of today's spirals. 

The epoch of the emergence of bars thus traces the 
epoch at which modern thin disks are established and 
start to dominate the stellar mass distribution in present- 
day spirals. This is primarily the interpretation of the 
bar fraction in our simulation sample, but if the red- 
shift evolution of the bar fraction agrees with observa- 
tions (which is the case up to z ~ 0.8, and needs to 
be probed at higher redshift) then it could suggest that 
the same interpretation applies to observed bars, namely 
that their emergence traces the typical epoch at which 
spiral galaxies establish their modern disk-dominated 
structure at z = 0.8 — 1 (for present-day stellar masses 
in the 10 10-11 M range). We will further discuss this 
hypothesis of a two-phase formation history traced by 
the emergence of bars, and illustrate with representative 
examples, in Section [7] 

5. BAR LIFETIME 

As shown in Section [3.11 bars start to steadily appear 
at redshift z ~ 0.8 — 1. A closer examination of in- 
dividual bars formed at this redshift reveals that their 
strength often remains roughly constant down to redshift 
z ~ 0. Is this because the conditions are favorable for the 
strength being intrinsically constant, or would those bars 
tend to weaken or even be destroyed in complete isola- 
tion but are maintained by e xternal infall as proposed by 
iBournaud fc Combesl ((2002)? We here examine the in- 
fluence of the external accretion on the lifetime of a bar 
in realistic cos mological context (unlike the idealized ac- 
cretion used by IBournaud fc Combesl (|2002D ). We select 
four galaxies, hereafter labeled Galaxy 1 to Galaxy 4, 
with representative bar strength histories in the main 
sample, and we run new simulations without external 
accretion once their bar has formed. 

Figure QT] shows the two first galaxies (Galaxy 1 and 
2), for which the redshift evolution of the bar strength 
is very similar with and without accretion. The bar is 
formed at z ~ 1.3 and z ~ 1 in these two galaxies, respec- 
tively. In both galaxies the bar strength is roughly con- 
stant down to z ~ with and without accretion. These 
bars are intrinsically long-lived and their evolution is not 
significantly influenced by late mass accretion. 

FigureQ2]shows the two other galaxies, for which exter- 
nal accretion has some influence on the bar strength evo- 
lution. In the fiducial simulation, Galaxy 3 (left panel) 
hosts a weak bar [S < 0.2 at all redshifts), that dissolves 
and re-forms several times. Without external accretion, 
the bar evolution is different, but the overall trend is sim- 
ilar, with several episodes of formation and dissolution. 
In this galaxy the bar remains weak in all cases. 

Galaxy 4 hosts a strong bar in the fiducial simulation, 



with S > 0.3 from redshift z ~ 0.8 down to z ~ 0. 
When we shut down external accretion, the bar stays of 
a similar strength down to z ~ 0.2, and then weakens 
significantly. In this case, external accretion is necessary 
to maintain a strong bar down to z = 0. 

Hence, bars that are formed after z ~ 0.8 — 1 are 
mostly long-lived, often without requiring external accre- 
tion, which is probably explained by the fact that these 
bars form once the gas fraction is relatively low - bar dis- 
solution is easier to reach with high gas fractions. There 
are however cases where the bar evolution is influenced 
by external mass infall and late cosmological accretion 
below redshift one, including cases where the presence of 
a strong bar at z = is only achieved through external 
accretion. 

At the opposite, the (relatively rare) bars that are 
formed early (1 < z < 2) have a shorter lifetimes. They 
can either dissolve spontaneously, as is the case for the 
first galaxy shown on Figure IT3l or be destroyed by merg- 
ers, as is the case for the two other galaxies shown on 
this figure. The shorter bar lifetime at z > 1 com- 
pared to z < 1 can be explained by the higher gas frac- 
tions at these early epochs, making easier to dissolve 
a bar by gravitational torquing and mass concentration 
(jBournaud et al.ll2005| ) . In these galaxies with early bars, 
the bar re-forms and eventually persist down to redshift 
zero, owing to late cosmological infall and build-up of a 
massive thin disk. 

Overall, it appears that external gas accretion is re- 
quired to maintain low-redshift bars only in a limited 
fraction of the sample, but is generally required to re- 
form early bars at z ~ 1, which tend to be intrinsically 
short-lived. 

6. THE ROLE OF BARS IN (PSEUDO)-BULGE GROWTH 

In this Section, we address the impact of the bar in 
the growth of the bulge, since bars have been proposed 
to play an important role in the formation of pseudo- 
bulge by trigger ing the gas inflows to the central re- 
gions of galaxy tiAthanassoulal 120051 : iHeller et all 120071 : 
iFisher et al.ll2009f ). 

He re the identification of th e bulge is made using GAL- 
FIT (jPeng et all 120021 . 1201(1 for galaxies at z = 0. Five 
galaxies (out of 33 analyzed so far) are removed from 
the study because their complex structure did not make 
it possible to achieve a satisfying and unique GALFIT 
decomposition. The majority of bulges have relatively 
low Sersic indices indicating tha t they are pseudo-b ulges 
rather than classical bulges (see iMartig et al.1l2012l for a 
detailed study of the simulated galaxies in our sample) . 

We consider all stars that are present in the bulge of a 
given galaxy at z = and we compute the averaged nor- 
malized distribution of the redshift formation of these 
stars. For galaxies that host a bar this distribution is 
computed with respect to the moment of the apparition 
of the bar which is thus interpreted as the average red- 
shift at which the bar forms. 

In Figure [H] we show the normalized age distribution 
of stars in the bulge for barred and unbarred galaxies 
(left and right panels, respectively), together with the 
theoretical ma ss accretion rate tha t scales with redshift 
as (1 + z) 2 - 25 (jNeistein et alj|2006f >. We note a slightly 
larger excess of stars in barred galaxies with respect to 
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Fig. 11. — Examples of two galaxies with similar redshift evolution of bar strength with and without accretion. In both cases the two 
galaxies host long-lived strong (5 > 0.3) bars that form at z > 0.8. The strength of bars does not evolve significantly. The external 
accretion hence does not seem to be needed to maintain a strong bar down to z = 0. Here and in all following figures, the null value for 
the bar strength is to be interpreted as no bar detection rather than exact value obtained by Eq. [2] 
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Fig. 12. — Examples of two galaxies with different redshift evolution of bar strength with and without accretion. Galaxy 3 (left panel) 
hosts a weak bar that forms at z > 1. Detailed redshift evolution of the bar strength in fiducial simulation is different than in the case 
without external accretion, but the overall behaviour is similar: bar undergoes several episodes of formation and dissolution. Galaxy 4 
(right panel) hosts a strong bar from z ~ 0.6 down to z ~ in fiducial simulation, but weakens significantly from z ~ 0.2 when external 
accretion is shut down. 

In our models, bars emerge in present-day spirals at 
z ~ 1, being rare and only weak at higher redshift, 
and rapidly becoming ubiquitous at lower redshift. We 
suggested from our simulations that this corresponds to 
the transition between an early "violent" phase down to 
z ~ 0.8 — 1 and a late "secular" phase at lower redshifts. 
Indeed, the morphology of the progenitors of today's spi- 
rals is rapidly evolving and uncorrclated with their final 
structure in the pr oposed "violent" phase at z > 1 (sec 
iMartig et al.ll2012l ). with mostly the thick disk and stel- 
lar spheroids forming in this phase, while after z < 1 the 
thin disk grows and the structural parameters such as 
the bulge and disk fractions become well correlated with 
their final values at z = 0. 

The influence of the dark matter halo can have 
additional effects on the evolution of bars (see e.g., 
Berentzen fc Shlosmanl !2006; M achado fc Athanassoulal 
20101 for the effects of the shape of halo and 



the infall rate as compared to that of unbarred galaxies. 

The qualitative comparison suggests that the presence 
of a bar has an impact on the growth of the bulge. We 
quantify this influence in terms of the ratio of the area of 
the histogram over the region between the redshift z = 
and the average redshift of bar formation to the area 
over redshift range of < z < 2. We find 24% for barred 
galaxies and 11% for galaxies without a bar. 

Even though it is difficult to quantify the effect of the 
bar on the bulge more precisely, both qualitative and 
quantitative study suggest that the bar plays a role in 
the growth of the bulge as cosmological accretion rate 
of baryons alone is not sufficient in explaining this over- 
abundance of star forming in the pseudo-bulge at late 
epochs. 



7. DISCUSSION: DOES THE EMERGENCE OF BARS AT 
Z ~ 1 TRACE A NEW PHASE IN SPIRAL GALAXY 
FORMATION? 



Weinberg fc Kate] [20071 : ISellwoodl [2001 for the effects 
involving more general halo properties). As the dark 
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Fig. 13. — Examples of three simulated galaxies with an early 
bar formation (1 < z < 2). Galaxy 5 (top panel) hosts a bar that 
is destroyed spontaneously rapidly after its formation at z ~ 2. A 
new bar reforms at low redshift (z ~ 0.2). The bars developed at 
z > 1 in Galaxies 6 and 7 (middle and bottom panels) are destroyed 
by mergers, but reform later. The early formed bars seem to be 
short-lived, although they can reform later. The high frequency 
at which the bar is destroyed and reformed between z > 1 and 
z < 2 in Galaxy 7 (bottom panel) is rather rare. Bars tend to 
undergo cycles of formation and destruction more commonly at 
higher redshifts, bars are also more rare at z > 1. 

matter halo evolves with redshift, it affects the evolu- 



tion of the entire galaxy, in particular the bar, which in 
turn influences the halo itself. Bar formation can thus 
be reinforced or delayed depending on exact halo proper- 
ties. Such effects should be resolved by our simulations, 
but the key epoch of bar formation appears to corre- 
spond mostly to the evolution of baryonic properties in 
our analysis. 

We show on Figure [15] three representative examples 
at z > 1: one, which is disk-dominated at z > 1 but 
where violent disk instabilities (including giant clumps) 
destroy the early disk into a thick disk and spheroid, 
another one, which is spheroid-dominated with a major 
merger at z ~ 2, and a third one, which is also spheroid 
dominated with several minor mergers at z ~ 1 — 2. This 
is illustrative of the "violent phase" at z > 1 . Indeed, it 
is known from other work that high-redshift disks have 
high gas fractions and are violently unstable with gi- 
ant clumps and tr ansient features, but do not frequently 
devel op bars (e.g.. iBournaud et alj|2007l : iCeverino et al.l 
I2010T ) and this instability destroys any thin disk that 
would have started to grow, while ma jor mergers that 
can r eform some disk components (e.g.. IRobertson et al.l 
120061) but mostly conve rt disks into spheroid s even in 
high-redshift conditions (jBournaud et al. 20TT[). 

Hence, in the two examples on Figure [T5] and in the 
majority of our sample, no massive thin disk can stabi- 
lize before redshift one and develop a substantial bar. On 
the other hand, mergers with mass ratios larger than 5:1 
arc almost absent from our sample after z — 1 and diffuse 
gas accretion occurs at much lower rate, and the rate of 
stellar bulge growth also drops. Thus a massive thin disk 
can form and start to dominate the mass distribution, as 
probed by the formation of spiral arms between 2 = 1 
and z = 0.5 in the face-on images shown on Figure [TSI 
This thin disk grows secularly down to z = and gen- 
erally gets barred by z « 0.5 - strongly in the two first 
cases, weakly in the third one. 

These three cases illustrate the transition between an 
early violent phase with frequent mergers and disk in- 
stabilities, and a late secular phase dominated by slower 
mass infall, and the fact that this transition occurs when 
the massive thin disk forms and is traced by the emer- 
gence of bars. Overall, the epoch of bar formation in 
our simulations probes the epoch at which spiral galax- 
ies have formed the bulk of their disk, stellar halo and 
thick disk, and start to be dominated (in stellar mass) 
by their final thin disk with only slow (secular) evolution 
down to z = 0. 

In observations, the bar fraction also decreases with 
increasing redshift, and although it could be probed in 
detail only up to z ~ 0.8 it is in close agreement with 
our models so far. In the following we discuss whether 
the emergence of bars at z ~ 0.8 — 1, if confirmed obscr- 
vationally, could also correspond to the epoch at which 
spiral galaxies acquire their final morphology and start 
being dominated by their thin rotating stellar disk. We 
are thus led to wonder whether a transition between a 
violent phase of galactic assembly, with rapid episodes 
of merging and violent instabilities building a thick disk 
and spheroids, could be followed by a calm secular phase 
of thin disk growth and evolution, with the typical tran- 
sition at z ~ 0.8 — 1. 

Chemical properties of disk galaxies could be con- 
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Fig. 14. — Normalized bulge stellar age distribution. The two panels show the age distribution for stars in the bulge at z = for different 
types of galaxies (left panel: barred galaxies, right panel: unbarred galaxies). The fitted theoretical total mass accretion rate that scales as 
(1 + z) 2,25 is shown too (the red line). The black dashed line represents the average time at which the bar is formed. All stars located to 
the left of this line are stars younger than the bar. We show this average bar formation rcdshift also for unbarred galaxies for comparison. 

itors of today's spiral galaxies in the 10 10-11 Mq mass 
range. This survey furthermore suggests that many of 
these z ~ 0.6 — 0.8 galaxies have formed their disks only 
recently after undergoing vio lent events such as major 
mergers (|Hammer et al.ll2009f ) and that they will undergo 
only slow evolution of their global proper ties, such as 
their Tully-Fisher relation, down to z = (|Puech et al.l 
[201?]) . 

Morphological studies are also consistent with the 
emergence of modern thin disks around redshift z ~ 
1. The Hubble Ultra Deep Field sample studied by 
lElmegreen et all (|2007al I2009D at z > 1 is dominated 
by irregular morphologies corresponding to major merger 
and interactions and "clumpy" unstable disks, which are 
typically forming thick disks, bulges and stellar haloes 
(|Bournaud et al.ll2009f ). not thin spiral disks. Some cold 
spiral disks are found in this sample but their fraction 
is quite low before 2 = 1. At z ~ 0.7, the situation 
is largely different as the fraction of clumpy irregular 
disks drops steadily and stable spiral disks rapidly be- 
come more numerous (still for progenitors of present- 
day IP 10 " 11 Mf) galaxies). Fo r somewhat more massive 
galaxies, iSargent et al.l (|2007l ) found that a large frac- 
tion of massive disks are in place around redshift one, 
but substantially fewer than at lower redshift. 

Hence, existing observations are consistent with our 
suggestion that the emergence of galactic bars at z~ 
0.8 — 1 traces the transition between an early "violent" 
phase during which stars that belong to the modern thick 
disk, bulge and halo arc formed in systems that do not 
have a permanent disk-dominated structure, and a late 
"secular" phase of thin disk growth and evolution with 
ubiquitous bars and limited pseudo-bulge gr owth. In this 
pictu re, the downsizing of bar formation (|Sheth et al.l 
120081 and Section 4.2) could correspond to the later ter- 
mination of the violent phase and later disk stabilization 
in lower-mass galaxies. This would be consistent with 
the fact that both merging activity and violent disk in- 
stabilities isJiojdd_jjej^ist^ i yjower redshift for lower mass 
galaxies (|Bournaud et al.l 120121 and references therein). 
It is also possible that bars grow more rapidly in more 



sistent with such two-phase assembly of their compo- 
nents. Thick disks are ubiquitous around spira l galaxies 
(jDalcanton fc Bernstemll2002t ISeth et al.ll2005D . without 
star for mation and no or almost no young stellar pop - 
ulations (|Yoachim fc Dalcantor]|2006l: ilbata et alj|2009f) . 
The Milky Way thick disk contains no or a l most no star 
youn ger than 8 Gyr (|Gilmore et al.|[l985t iReddv et al.l 
120061 ). This is consistent with the occurrence of events 
disrupting any cold rotating disk into a thick disk only 
at redshifts higher than z ~ 1, while the thin disk 
wouldn't have been significantly disrupted/thickened af- 
ter this epoch. The enhanced a el ement abundances in 
the thick disk and central bulge (|Lecureur et al.l [20071 : 
iZoccali et al.ll2007j) suggest that their star formation oc- 
curred mostly in brief events (not longer than a few 
10 8 yr), and that the formation of stars that belong to 
the present-day thin disk occurred at late r epochs with 
longer timescales (see also iChiappinil 120091 ) . These brief 
events could have been mergers, as well as violent in- 
stabilities (giant clumps) in rapidly-ac creting galaxies, 
that grow both a bulge and a thick disk (jBournaud et al.l 
120091) over short timescales, while a thin disk component 
can form and remain stable after these phases with lower 
gas infall rates. 

Direct searches for disks, through gas kinematics 
and/or optical and near-infrared spectroscopy, are also 
consistent with z ~ 1 as the typical redshift for the 
emergence of modern thin spiral disks, along with their 
bars. Near-infrared spectroscopy of a z ~ 3 sample of 
star- forming galaxies (AMA ZE and LSD. lMaiolino et al.1 
[20101 IGnerucci et al.|[201lT) . with masses typical for the 
progenitors of Milky Way-like galaxies and present-day 
spirals, reveals only a minor fraction of rotating disks, 
the majority of these galaxies being mergers or irregu- 
lar systems dominated by high velocity dispersions. A 
similar survey of star-forming galaxies in a compara- 
ble mass range at a median redshift z ~ 1.2 (MASSIV, 
lEpinat et all 120 121 ) finds about 40-50% of rotating disk 
galaxies. At z ~ 0.6 - 0.8, the IM AGES survey find s a 
majority of rotating disk galaxies ()Yang et al.ll2008l) in 
a sample that still covers masses typical for the progen- 
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Fig. 15. — Examples of morphological evolution from 2 = 2 down to z = for three simulated galaxies. Stellar density maps (face-on 
50 X 50 kpc 2 projections) are shown for z = 2, 1, 0.5 et 0. Galaxies evolve rapidly during the "violent" phase at z > 1 when they 
frequently undergo phases of violent disk instabilities (top panels) and major mergers (middle panels). The morphology at z = and 
z > 1 are uncorrclatcd with early spheroids being possible progenitors of today's spirals (bottom panels). Once galaxies enter the late 
"secular" phase at z < 1, their structural parameters become more tightly correlated with the final disk morphology. The color coding of 
the projected density maps is the same as in Figure [T] 



massive systems once th eir cold, thin disk is stabilized 
(|Elmegreen et al.ll2007bO . The observed evolution of the 
bar fraction so far is consistent with our model, but fur- 
ther confirmation could be obtained by confirming the 
drop in the bar fraction at z ~ 1 and above, with almost 
only weak bars (strength < 0.2) being present at z > 1 
for the mass range studied here. 

8. SUMMARY 

We studied a sample of cosmological zoom-in simula- 
tions of 33 Milky Way-mass galaxies in field and loose 
group environments from 2 = 2 down to z = 0. The 
method used to determine the presence of a bar is based 
on the decomposition of the stellar surface density pro- 
files into Fourier components. We also analyzed the 
disk/spheroid structure of the modeled galaxies using the 
Sersic index of the surface density profile. 

Our main results are as follows: 

i) The total bar fraction among spiral galaxies declines 
with increasing redshift. It drops from almost 90% 
at z = to about 50% at z ~ 1 and to almost zero 



at z ~ 2. The fraction of observable and strong bars 
declines from about 70% at z = to 10 — 20% at 
z ~ 1 and to zero at z = 2. This result holds for 
galaxies with mass range of 2 x 10 9 — 8 x 1O 1O M0 at 
z ~ 1 and 4 x 10 5 - 1 x 1O 1O M at z ~ 2, i.e. typ- 
ical progenitors of Milky Way-like spirals. However, 
the bar fraction could remain higher at z ~ 2 for 
more massive galaxies, if the downsizing of bar for- 
mation observed in our sample still holds for higher 
masses /redshifts. 

The epoch of bar formation traces the epoch of the 
emergence of final disk of spiral galaxies. This cor- 
responds to the termination of an early "violent" 
phase at z > 1, characterized by frequent mergers, 
violent disk instabilities and rapidly evolving struc- 
ture, forming thick disks, bulges and stellar halos. 
It is followed by a "secular" phase at z < 0.8, domi- 
nated by the slower growth and evolution of modern 
thin disks and limited bulge growth at late times. 
The z = 0.8 — 1.0 transition epoch is for the mass of 



14 



typical Milky Way progenitors, and tends to move 
to higher rcdshift for more massive systems. 

iii) We find that there is only a minor contribution of 
bars in the late growth (z < 1) of (pseudo-) bulges 
in spiral galaxies. This late growth is dominated 
by continued cosmic infall and minor mergers rather 
than by bars. 

iv) Finally, early bars (formed at z > 1) are often short 
lived and may reform several times. Bars formed be- 
low z ~ 1 are found to persist down to z = 0, some of 
them being intrinsically short-lived but maintained 
by late cosmological gas infall. 

If confirmed obscrvationally, the scarcity of significant 
bars at z > 1 would indicate, according to our mod- 
els, that present day spirals and Milky Way-like galaxies 
have formed and stabilized their modern thin spiral disk 
only relatively late in their growth history, typically at 
z ~ 0.8 — 1. At earlier times, they would be mostly form- 



ing their spheroidal components (bulges, halos) and thick 
disk, under the effect of both hierarchical merging and vi- 
olent instabilities in rapidly-accreting systems. The con- 
tinuation of this violent phase with mergers, rapid cold 
gas accretion and disk instabiliti es down to lower red- 
shift in lower mass galaxies (e.g.. iBournaud et alJl20i2T ) 
could then explain a "downsizing" -like behaviour for bar 
formation. 
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